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Abstract 

A prescription is presented for real-time finite-temperature perturbation theory in co- 
variant gauges, in which only the two physical degrees of freedom of the gauge-field 
propagator acquire thermal parts. The propagators for the unphysical degrees of free- 
dom of the gauge field, and for the Faddeev-Popov ghost field, are independent of 
temperature. This prescription is applied to the calculation of the one-loop gluon self- 
energy and the two-loop interaction pressure, and is found to be simpler to use than 
the conventional one. 
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1. INTRODUCTION 



In finite-temperature field theory, one calculates a grand partition function 

i 

and thermal averages 

(Q) = Z-i^(i|e-'5^g|i). (1.16) 

i 

In both cases, the sum is over a complete orthonormal set of physical states. In scalar field theory, these 
states span the whole Hilbcrt space, and so one may replace the sums with traces of the operators. 
But in the case of a gauge theory the Hilbert space contains also unphysical states and so things are 
more complicated. 

In this paper, we reconsider what to do about this for the case of covariant gauges. The traditional 
solutionis] is, in effect, to sum over all states of the gauge particles, and then cancel the unphysical 
contributions with ghosts. We shall argue that, for calculational purposes, it may be simpler to include 
in the sums (1.1) only the physical states from the start, and not to introduce unphysical terms that 
have to be cancelled. 

The finite-temperature propagator consists of a vacuum piece, plus a "thermal" part. In real-time 
perturbation theory, which we consider initially, it is a 2 x 2 matrix^] . Its vacuum piece is 

,•r,M<^^,^ _ ( {0\TA'^ix)A-iO)\0) (0|^^(-ia, O)A'^(x)|0) \ 

' W - l^(0|^M(a;0 _i(^^x)yl-(0)|0) (0|f.4'^(x)y4''(0)|0) J ^ ' 

with T the anti-time-ordering operator, and < a < /3. There is a similar matrix ghost-field propa- 
gator. The elements of the thermal part of the matrix take account of all the possible states of the 
heat bath and involve the Bosc distribution. In the traditional approach to covariant gauges, all the 
components of the gauge-field propagator have a thermal part, as also does the ghost propagator. 
We shall argue that it is simpler to give only the propagators for the physical degrees of freedom a 
thermal part, so that the unphysical components of the gauge field, and also the ghost field, have only 
vacuum propagators. In the next section, we give a derivation of this from the Gupta-Bleuler method 
of quantising the gauge field, though in a sense the result is obvious and can be written down without 
any long derivation. 

This is because the thermal part of a propagator contains a (5-function that puts the corresponding 
particle on shell. The thermal part represents the absorption of a particle from the heat bath or the 
emission of one into it. That is, the thermal part is directly associated with the particles in the physical 
states I i ) summed in (1.1). Consequently only the physical degrees of freedom acquire thermal parts 
for their propagation. 

The thermal part of the propagator for the gauge field therefore has the tensor structure 

Mk)er{k) (1.30) 



A=l 

where ex{k) are polarisation vectors and k'^ = 0. There is some freedom in the choice of these, but 
since the heat bath breaks the Lorentz symmetry already, the obvious choice is the one that causes no 
additional Lorentz-symmetry breaJcing. So we choose polarisation vectors orthogonal to the 4-velocity 
u of the heat bath, and then (1.3a) is 

u.k (u.ky 

In the frame where the heat bath is at rest this vanishes if either or i/ is zero, and its other elements 
are 

r.(fc) = j^.___. (1.3c) 



To summarise what we have said, part of the finite-temperature propagator of the gauge field is the 
vacuum piece, which in Feynman gauge reads 



l,-27ri<5(+)(fc2) -l/ij^-ie) )■ ^^'^''^ 

We have not written exphcitly the colour factor d""^ . We have made the particular choice a = in 
(1.2) so that the vacuum part of the propagator is temperature- independent; for other choices of a it 
depends on /3, but in a simple way. To (1.4a) must be added the thermal piece of the propagator: 

-iT>^%k)2T:5{k^)n{\ko\)(^^ , (1.46) 

where n{x) = l/(e^^ — 1). 

There is also the ghost propagator; this has only the vacuum part, which is just equal to the matrix 
that appears within (1.4a). Wc derive the propagator (1.4) in section 2 and generalise it to other 
gauges, both covariant and noncovariant. 

In order to show that this prescription can be simpler than the conventional one, we use it in section 3 
to calculate the gluon self-energy and the plasma pressure to lowest non-trivial order in the coupling. 
Section 4 is a discussion of our results. In an Appendix we introduce some alternative versions of the 
formalism in which the 2x2 propagator matrix is diagonal and temperature-independent, and the 
dependence on the temperature is instead transferred to the interaction vertices. 



2. GUPTA-BLEULER THEORY AT FINITE TEMPERATURE 

To evaluate the grand partition function (1.1a) or a thermal average (1.1b), we need to express the 
Hamiltonian in operator form and to identify the physical states. For this, we use standard Gupta- 
Bleulcr thcoryt"^]. Wc go through the details only for the case of the Feynman gauge. 

Start with the standard QCD Lagrangian — j-F^, with the addition of a gauge-fixing term but initially 
no ghost part: 

>C = -iF^^F^.-|(a.A)2 (2.1) 
The canonical momenta conjugate to the space components of the gauge field are 

tt' = (2.2a) 

while that conjugate to A'^ is 

7r° = -d.A . (2.26) 

The Hamiltonian density is 

H = 7r°i° + w'A' -C = Ho + H^^^. (2.3) 

Here, Hq is the part of H, expressed as a function of the fields and the canonical momenta with no 
explicit time derivatives of the fields, that survives if the coupling g is set equal to zero, and 

i?™^(7r,^) = -gn'.A'^ aA + ^{F'^F'^Y^^. (2.4) 

To derive perturbation theory, one introduces interaction-picture operators Qi{t) = A{t)Q{t)A~^(t), 
where A(t) = (=^tioi{t-tQ)^-zH{t-to) ^ These operators evolve with time as free fields: they obey Hamil- 
ton's equations of motion with Hamiltonian Hqi and the relation between the interaction-picture 
canonical momenta and the fields is as in (2.2) but with the coupling g set equal to 0: 

ni = Aj - d'A^i, TT? = -d.Aj. (2.2c) 



The interaction-picture field equations are just the zero-mass Klein-Gordon equations. Their solutions 
are 

(^^e-^'-^«''W + hc (2.5) 

where = |k|. The canonical equal-time commutation relations [Aj(i,x),7rj(i,y)] = iS^'^ 5^^\:k — y) 
etc give 

[ai'{k),a''^k')] = -5'^^2|k|(27r)3(5(3)(k-k'), 

K(fc),a'^(fc')] = [a^hk),a''hk')] = 0. (2.6) 

We now have to identify the physical states. For finite values of Iq, the time at which the interaction 
picture coincides with the Heisenberg picture, the constraint that they satisfy is non-trivialW . How- 
ever, to derive zero-temperature perturbation theory, or finite-temperature field theory in the real-time 
formalism, one may switch off the gauge coupling adiabatically and work with the limit to — > — oo. 
Then the interaction-picture states are just the in-states and the condition that picks out the physical 
states is the same as in abelian gauge theory. Its most general form is that matrix elements of d.Aj 
vanish between physical states, but we may impose the stronger condition!"^] 

a°{k) I phys, in) = = k.a{k) \ phys, in ). (2.7) 

In particular, the vacuum is required to satisfy these conditions. Together, (2.5), (2.6) and (2.7) make 
the vacuum expectation values of T-products of the fields take the usual Feynman-gauge-propagator 
forms. 

The time-development operator for the interaction-picture states away from t = to is A(t). By con- 
verting the differential equation it satisfies into an integral equation and iterating, one finds as usual 
that the S-matrix is a time-ordered exponential of an integral of H^^"^ {Aj ,ttj). if™""" is given in (2.4). 
At this stage we may replace the ttj with their expressions (2.2c) in terms of Aj and Aj. This gives 

^INT ^ _£lNT(^^^^^^ ^ ^g\A'jAA})\ (2.8) 

The last term here compensates for the presence of gA\.A^ A A^j in the interaction Lagrangiant^l . 
The Hamiltonian formalism of perturbation theory that we have outlined requires that, when 
the A\ field propagates between two such neighbouring vertices one should use the propagator 
{0\TA\{xi)A^j{x2)\0)- However, the usual Feynman rules use instead the double time derivative of 
the ^l/'field propagator, where the two time differentiations are applied also to the ^-functions that 
appear in the definition of the T-product. The last term in (2.8) is then simply omitted . 

From the differential equation satisfied by A(t) it is straightforward to show that the 5-matrix is 
unitary. However, this does not imply that probability is conserved, because the completeness relation 
includes the unphysical states. In order to achieve probability conservation one must add to the 
Lagrangian (2.1) a ghost part. This then guarantees that physical initial states do not scatter into 
unphysical final states. One may show this directly in the operator formalism^, though we shall be 
content here to accept the usual Faddeev-Popov path integral arguments as justification. 

We now extend the perturbation theory to nonzero temperature. The same operator A enters again, 
now with complex argument. Simple algebra gives I*^] 

Z = ^{i |e-/5^-A(to - iP)A-\to)\ i ) 

i 

{Qit)) = Z-'^{i |e-^^°^A(io - i/3)A-\t)Qj{t)A{t)A-\to)\ i ) (2.9) 

i 

where to must be the time at which the interaction picture coincides with the Heisenberg picture and 
the sum is still over all physical states, which we now choose to be interaction-picture m-states. 



In order now to derive the Feynman rules, it is necessary first to check that Wick's theorem may be 
used as usually in thermal field theory. For this, note that the physical-gluon operators commute with 
the unphysical gluon and ghost fields. Also, Hqi is a sum of commuting parts containing respectively 
only physical and unphysical operators, Hqi = H^^^^ + H^^^^^^ , and H™^^^^ has zero eigenvalue 
in the physical states. Hence any unphysical operators in the product of operators in (2.9) may be 
factorised out: only the vacuum expectation value of their product appears, and the normal zero- 
temperature Wick's theorem applies to this. For the remaining factor, we may replace e~^^"' by 
Q-I^^Qi ; the normal finite-temperature Wick's theorem applies to this factor. 

So the real-time finite-temperature perturbation theory is much as usualP], except that for the un- 
physical fields the propagators are just vacuum expectation values (1.2), both for the unphysical 
components of A and for the ghost. 

The gauge-field propagator in Feynman gauge therefore reads 



D^, - g„u g-afco0(jto) (Af - A^) -A|, 

1 / 1 e'^''° \ 



F) 

(2.10) 



with < (T < /3, and Aj? = l/(/c^ + ie). The choice a = (3/2 is usually made on grounds that this 
makes the real-time propagator a symmetric 2x2 matrix. Here this is the case only for the transverse 
projection T^^^Y) p„T"'^ . A more natural choice is o" = 0, for it renders the vacuum piece independent 
of (3. In the Appendix, we discuss various options for the diagonalisation of (2.10). 

Although the derivation was done in Feynman gaTigc, it seems obvious how the resulting propagator 
will read in other gauges. In a general linear gauge with quadratic gauge breaking term ■^{A^^ ffj^f^A"), 
where is a 4- vector which is either constant or constructed from derivatives d/dx, the vacuum piece 
generalises by 

g,. - G,.{k) = g,. - Mji±Mj± + (/2 _ Cfe^)^. (2.11) 

Here / = f{ik), if the gauge breaking term contains derivatives. With = ik^, this reproduces 
the propagator in covariant gauges, while = (0, ik) corresponds to the Coulomb gauge. The ghost 
propagator is determined by A~j^^^g^ = /./c; its 2 x 2 matrix structure is analogous to the vacuum 
part of the gauge propagator. In most gauges, other than Feynman or Coulomb, an infinitesimal 
nonhermitian piece has to be included in G^^ in order to give meaning to the denominators in (2.11). 
In this case, the vacuum part of D^;^ becomes 



e-^oe{ko) [ApG,, a>g;,) -a;.g;. 



(2.12) 



Although, strictly speaking, we have given a derivation only for the Feynman gauge, the results we 
have quoted agree with those derived by James for the = gauge!'']. For general covariant gauges, 
they agree with those in the literature!^] if one sets the temperature to zero in the non-T'*'' part of 
the propagator. 



3. CALCULATIONS 

3.1. One- loop gluon self-energy 

In a high-temperature expansion, the leading-temperature contributions arise from one-loop diagrams, 
and they determine the physics at the scale gT, where g is the coupling constant. In particular, 
the one-loop self-energy determines the spectrum of quasiparticles in the plasma, and contains the 
information on screening and Landau damping. It has been shown by explicit calculations P] and 
also by an analysis of the relevant Ward identitiest-*-^] that these "hard thermal loops", as they have 



been called, are independent of the gauge conditions needed to define the gluon propagator. With 
the usual calculational procedure, this gauge independence arises in a rather nontrivial manner, and 
in particular in covariant gauges the contributions from the thermalised Faddeev-Popov ghosts are 
decisive. With our prescription, the latter are absent and it turns out to be much simpler to establish 
the gauge independence. 

With the usual prescription for covariant gauges, the Feynman diagrams containing the hard-thermal- 
loop contributions to the gluon self-energy are those given in figure 1, where a slashed line denotes 
the thermal part of a propagator, and a bare line the vacuum part. With both vertices of type 1, 
these diagrams give Re ILn which coincides with the real part of the proper self-energy correction 11 
for timelike momenta. 

In our formalism, we do not have contributions from the ghost diagrams Ic, Id, although in diagram la 
the gauge-dependent vacuum gluon propagator does contribute. We shall show that the result is 
nevertheless gauge independent. 

With our real-time Feynman rules, the diagram lb is identical to the corresponding one in Coulomb 
gauge. Since in Coulomb gauge there are also no contributions from the ghost diagrams, diagram la 
has to be shown to coincide with the corresponding one in the Coulomb gauge. Evaluating diagram la 
in Coulomb gauge, one finds for its leading-temperature part in SU(N) gauge theory 



I) 



-iD-2) 



Ak^ky Ik^Qi, iQukn 



A-2 



(k-qy 



;T^u{k) 



(3.1) 

where the principal value of the pole of the integrand is understood. Apart from the piece containing 
the Feynman denominator, only the terms with highest degree in k contribute. For this reason, one 
may further simplify k^/^k. — q)^ ^ 1 in the last term of (3.1). 

On the other hand, with the Feynman-gauge propagator (2.10) one obtains 



^iii<l)=9'Nj^^0^d{k^)n{\ko 



I) 



[k-qf 



{k-qf 



(3.2) 

Because of the presence of the delta- function d{k^), the numerator in the last term may be rewritten 
as {k + q)^ —{k — q)^ + 2q^, and the dropped because it does not contribute to the hard thermal 
part. Thus (3.2) gives the same leading-temperature terms as (3.1) does. 

Adding the contribution from diagram lb gives the well-known result 



(27r)(^-i) 



S{k^)n{\ko\ 



9nu 



Akfj^ki, Ikpi^qi, Iq^kn 



(3.3) 



In gauges other than the Feynman gauge, the vacuum part of the gluon propagator has to be changed 
according to (2.11). Diagram lb is not affected by this, but diagram la receives an additional contri- 
bution 

d^'k 2//^ ^ {^^{k-qf-f){k-qf- 



g^N 



(27r)(-D- 



^6{e)n{\ko\)TP%k) 



+ 



[/•(fc-9)P 



(3.4) 



X (fe - q)"Vcij,aik - q, q, -k)Vp^i3{k, -q, q-k), 



where in general = /^(A; — q). V is the 3-gluon vertex (without the SU(N) colour factor), which 
obeys 



{k - qYVa^a{k - q,q, -k) = gafj,{k -q )- k„k^ + q^qf,. 



(3.5) 



When this is inserted into (3.4), only the two terms quadratic in the loop momentum k contribute to 
the leading-temperature part, but both get cancelled, the k"^ by the (5-function, and the kak^i by the 
transverse projector Tpcr{k). Therefore the result (3.3) is completely independent of the gauge fixing. 

For completeness, we also give the result, in our formalism, for the one-loop contribution to the 
damping of long- wavelength transverse gluons. In the conventional formalism, the imaginary part 



of the gluon self-energy at one-loop order is given by the diagrams of figure 2 and is strongly gauge 
dependent. In covariant gauges, the result comes with an unphysical sign, which was even taken to be 
suggestive of a fundamental instability of the perturbative thermal ground state (see reference 11 for a 
review). In our formalism, only diagram 2a contributes, and the result is obviously gauge independent; 
it reads 

lm{T^,U>^^){ko,0) = -^Ng^T\ (3.6) 

By virtue of Im(T^j,n^^) = tanh(fco/3/2)Im(T^j,n^^) this gives a contribution to the damping rate of 
long-wavelength gluons 

2 n~i AT 

one-loop _ 9 J-^^ ^\ 

' 247r ■ ^ ' ' 

This result coincides with the Coulomb-gauge one of the conventional formalismt^^l , where also only 
transverse gluons contribute to the imaginary part. It is well-known by now 

[12] 

that one-loop results 

for the gluonic damping rate are incomplete, and the gauge independence of our result (3.7) should 
not lead one astray. Still, we find it gratifying that in our formalism the unphysical modes are not 
able to contribute to (3.7), whereas normally in the bare one-loop calculation they even give rise to a 
negative sign. 



3.2. Two-loop gluon-interaction pressure 

The thermodynamic pressure of the gluon plasma, P = 1/{PV) In Z, is a physical quantity and there- 
fore gauge independent. In the following we shall calculate it up to two- loop order in order to perform 
a further check on our Feynman rules. The diagrammatic rules for calculating the partition function 
in the real-time formalism have been described in reference 13, to which we refer for more details. 

At one-loop order, our formalism leads to exactly the same expression for P as the conventional one 
does in Coulomb gauge. Differences occur starting at two-loop order. The diagrams to be considered 
are given in figure 3, where thermal contributions have to be inserted in all possible ways, keeping at 
least one of the vertices appearing therein of type 1 

[13] 

Diagram 3b potentially gives rise to integrals involving three powers of the Bose distribution function, 
but these we find to vanish after performing the momentum algebra because of the three delta functions 
associated with the distribution functions, exactly as in the conventional formalism. 

Next consider the diagrams where two lines are thermal. Diagram 3a then identically reproduces the 
result of conventional Coulomb gauge, 

P^^ = - iAr(iV^ - 1)/ / '^^n{\qo\)6{q'M\ko\)S{k^)2 (3 - z^) 

J (2^) (3.8) 

= - ^^(^' - ^)9'T\ 

where we introduced z = k.q/|k| |q|. Diagram 3b with two thermal and one vacuum gluon propagator is 
different, though. The vacuum piece of the gluon propagator contributes potential gauge dependences. 
However, because of (2.11), all gauge-dependent terms have one 4-momentum contracted with a 3- 
gluon vertex, and give rise to integrals of the form 

j d\d''kn{\qo\)6{q^)n{\ko\)S{e)T^'-{k)T^f{q){k - qrV^aa{-K k-q,q)---. (3.9) 

Because of 

{k - q)"Vf^aa{-k, k-q,q) = g^^iq'^ - k'^) - qaQti + Kk^, (3.10) 



it is easy to see that all the integrals of the form (3.9) vanish — (f and k'^ in (3.10) disappear because 
the (5-functions in (3.9), and the other terms also owing to the presence of the transverse projectors 



T'^'^{k)T'^P{q). Thus the contributions from figure 3b which involve two powers of distribution func- 
tions are completely gauge independent. There are actually two such diagrams. The one with both 
3-vertices of type 1 gives 



(1 + z^){k + qf + 8(1 - z2)(k2 + q2) ^^•^^'^^ 



where the first term can be simplified according to {k + q^ /{k — q)"^ —1 because of the delta 
functions, and the last term vanishes after performing the integrals, making this contribution coincide 
with the one of conventional Coulomb gauge. This yields 

Because the vacuum part of the gluon propagator (2.10) does have a (12)-component, there is a second 
diagram of the form 3b, where one vertex is of type 1, and the other of type 2. We find after carrying 
out the momentum algebra that the latter vanishes because of the presence of three (5-functions, 
S{q'')Sik')S{{k-qn 

Finally, there are the diagrams with only one gluon propagator thermal. Those diagrams, however, 
contain a zero-temperature one-loop subgraph with its external lines put on-shell by the thermal 
propagator, and so are removed by renormalisation. The additional diagrams containing both a type- 
1 and a type-2 vertex do not complicate this picture, for they turn out to vanish identically. 

Discarding the purely zero-temperature contributions, we have therefore that the gluon-interaction 
pressure at two-loop order is given by the sum of (3.8) and (3.11), 

P2 = -^N{N^-l)g^T\ (3.12) 

reproducing the standard result I ^'^"''^^1. Notice that the ghost diagram 3c did not contribute to (3.12), 
although we worked in a general linear gauge. 



4. CONCLUSION 

We have demonstrated that in the real-time version of finite-temperature perturbation theory it is 
quite natural to include a thermal part in the gauge propagator only for the two physical degrees of 
freedom of the gluons, leaving the unphysical gauge-field components and the Faddeev-Popov ghosts 
unheated. This allows us to combine the simplicity of covariant gauges for the zero-temperature 
part with the advantages of the noncovariant ghost-free gauges for the thermal contributions. For 
the thermal contributions we single out the rest frame of the heat bath and do not cause additional 
violation of Lorentz symmetry. Our derivation has been carried through for the Feynman gauge. 
The generalization to arbitrary gauges is immediate, though there may be complications from the 
prescriptions necessary to define the poles at f.k = in (2.11). 

We have successfully tested these new Feynman rules in one- and two-loop calculations, and in par- 
ticular have been able to demonstrate gauge independence of the high-temperature part of the gluon 
self-energy and for the 2-loop gluon-interaction pressure in a remarkably simple (albeit nontrivial) 
and general manner. In the conventional formalism, gauge independence of the hard-thermal gluon 
self-energy has been explicitly checked only in certain classes of gauges t^l, whereas the 2-loop pressure 
has been calculated so far only in Feynmant-'^'^'-'^^] and axialt-*^^'-*^^] gauges. 

We have worked throughout in the real-time formalism. This enables us to develop the perturbation 
theory in terms of in fields, with the great advantage that not only do the fields obey free-field 
equations, but also we may ignore the interaction in the condition that picks out the physical states. 



We have not been able to translate our prescriptions satisfactorily into an imaginary-time formalism. 
In order to set up such a formalism, one would need to choose an interaction picture that coincides 
with the Heisenberg picture at finite time to- The condition that picks out the physical interaction- 
picture states is then much more complicated!^]. Even if we were able to implement it, the fact that 
the unphysical-field propagators are temperature-independent and so not periodic in imaginary time 
would lead to extra complications!^] and make the formalism difficult to use. 

An interesting question, to which we hope to return in future work, is how the present approach 

could be extended to accommodate a resummation of finite-temperature perturbation theory, as is 
mandatory for exploring physics at the scale g^T in high-temperature quantum chromodynamicst^'^] 
(and, incidentially, for handling pinch singularities^"'^''']). An important difference between the bare and 
the resummed theory is that, because of the appearance of an additional collective mode at the scale 
gT, (1.3) no longer covers all of the physical modes in the plasma. Including the latter in a formalism 
analogous to the one presented here might prove to be of conceptual interest for the resummation 
program!-'^], and might perhaps allow an easier verification of its gauge independence. 

One of us (PVL ) is grateful to Peter van Nieuwenhuizen for patient discussions of the Gupta-Bleuler 
formalism; AR wishes to thank Patrick Aurenche and Randy Kobes for useful remarks. 

APPENDIX: Diagonalisation of the matrix propagator 

The real-time version of finite-temperature field theory can be reformulated by diagonalisation of the 
2x2 matrix propagators, which allows one to associate all thermal contributions with the vertices. 

A well-known possibility^' ^'^] is to diagonalise to a matrix constructed from the ordinary Feynman 
propagator Ap, 

In our formalism, the Feynman-gauge propagator may be written as 

n^"" = T^^ MtDMt - {g^"" + T^"") MqDMo {A.2a) 

with (when a = 0) 

and Mo the zero-temperature limit of Mt: 
The ghost propagator is MoDMq. 

Recently, a different scheme has been proposed by Aurenche and Becherrawyl-*-^] , which diagonalises 
to a matrix composed of retarded and advanced propagators, 

instead of Feynman propagators. For this, one has to introduce different matrices to be associated 
with incoming and outgoing lines [1^1. In our formalism for the Feynman gauge 



D^^ = T^""" UtDVt - {g^"" + T'^'') UoDVo 



{AAa) 



with 



The ghost propagator is UqDVo. The matrices U and V are removed from the propagators and 
instead associated with the in- and outgoing hnes of the vertices. (Because of 



""""Vi -(i + n(M)J' ^""li -e{k,), 

l + n{ko) n{ko)\ _fe{ko) -e{-ko)^ 



Vik)T^=V'{-k), n 



1 

1 



more symmetric Feynman rules result if a factor ri is combined with U and D, turning the latter 
off-diagonal[19] .) With (A.4), the analysis of reference 18 carries over to a large extent. However, in 
addition to retarded and advanced lines one has to distinguish between transverse and non-transverse 
ones, and in particular there is no simple causality principle when transverse and non-transverse lines 

come together at a vertex. For example, a 3-vertex connecting solely retarded or solely advanced lines 
does not vanish, unless the lines are either all transverse or all non-transverse. 

Such a reformulation in terms of retarded and advanced Green functions leads to well-defined causal 
properties, which are not manifest in the usual real-time formalism I^^l . Causal Green functions, 
which are relevant for example in linear response theory, are usually obtained in a direct manner in 

the imaginary-time formalism. In our prescription, however, this is not straightforward to implement, 
and so this diagonalization scheme appears to be particularly useful for our approach. 
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Figure captions 

1 Diagrams which contain the hard-thermal-loop contributions to the gluon self energy. Solid lines are 
gluons and broken lines Faddccv-Popov ghosts. The slashes denote thermal parts of the propagators; 
unslashed lines are vacuum propagators. In the formalism described in this paper, diagrams (b) and 
(c) do not occur. 

2 The imaginary part of the gluon self-energy, with the same notation as figure 1. In our formalism, (b) 
does not occur. 

3 Diagrams for the 2-loop gluon interaction pressure. The lines represent complete propagators (vacuum 
plus thermal parts). 



